Circular dichroism simulated spectra of chiral gold nanoclusters: A dipole 

approximation 
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Circular dichroism (CD) spectra of chiral bare and thiol-passivated gold nanoclusters have been 
calculated within the dipole approximation. The calculated CD spectra show features that allow us 
to distinguish between clusters with different indexes of chirality. The main factor responsible of the 
differences in the CD lineshapes is the distribution of interatomic distances that characterize the 
chiral cluster geometry. These results provide theoretical support for the quantification of chirality 
and its measurement, using the CD lineshapes of chiral metal nanoclusters. 
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Recent theoretical studies on the structural proper- 
ties of gold clusters have shown that the lowest-energy 
(most stable) isomers of bare Au28 an d AU55, and 
thiol-passivated Au28(SCH 3 )i 6 and Au 38 (SCH 3 )24 clus- 
ters correspond to chiral nanostructuresi*2. These results 
provided theoretical support for the existence of chiral- 
ity in metal clusters, suggested by the intense optical ac- 
tivity measured in the metal-based electronic transitions 
of size-separated glutathione-passivated gold nanoclus- 
ters in the size range of 20-40 atomsi 

Another theoretical prediction, based on quantify- 
ing chirality through the Hausdorff chirality measure 
(HCM)i2i4, indicates that the strong structural distor- 
tion in a gold cluster upon thiol passivation could induce 
chirality in an achiral bare cluster or increase the index 
of chirality in an already bare chiral clusteri^. To con- 
firm this prediction further experimental studies on chiral 
bare and passivated clusters are expected to be performed 
in the near future using, for example, circular dichroism 
(CD) spectroscopy^. 

To explore the capability of the CD technique to detect 
different indexes of chirality existing in bare and pas- 
sivated gold nanoclusters it would be useful to have a 
theoretical estimation of their CD lineshapes. Quantum- 
mechanical calculations of CD spectra for small and mod- 
erately sized organic molecules had been performed 6 , 
however, this type of approach is still inappropriate for 
passivated gold nanoclusters due to the huge computa- 
tional effort that would be involved 6 . In this work, we 
present results on simulated CD spectra of chiral gold 
nanoclusters obtained through a dipole approximation^. 
This study allows to gain insights on the variations of the 
CD lineshape due to different indexes of chirality exist- 
ing in gold clusters. It is expected that this information 
would be useful to motivate further experimental stud- 
ies that confirm not only the existence of chirality in gold 
nanoclusters, but also to correlate distinct features of the 
CD spectra with the different indexes of chiralitjiii^. 

In this work the nanoclusters of interest are composed 
by N gold atoms. We suppose that each one of these 
atoms is represented by a polarizable point dipole lo- 



cated at the position of the atom. We assume that the 
dipole located at r,-, i = 1, 2, . . . , N, is characterized by 
a polarizability cxi(u>), where ui denotes the angular fre- 
quency. The cluster is excited by a circularly polarized 
incident wave Ej nc = Eo(x±ty)e lk ' r ~ luJt , where t denotes 
time, k is the wave- vector, whose magnitude is given by 
k = uj/c = 2ir/\, c is the speed of light, and A is the 
wavelength of the incident light. Here x and y denotes 
the unitary vector along any of two orthogonal directions, 
and the sign + (— ) corresponds to right- (left-) circularly 
polarized light. 

Each dipole of the system is subject to a total electric 
field that can be divided into two contributions: (i) the 
incident radiation field, plus (ii) the radiation field re- 
sulting from all of the other induced dipoles. The sum of 
both fields is the so called local field and is given by 
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where Pi is the dipole moment of the atom located at r j , 
and is an off-diagonal interaction matrix with 3x3 
matrices as elements, such that 
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On the other 



hand, the induced dipole moment at each atom is given 
by Pi — a i ' E.^ioc, such that 3A-coupled complex linear 
equations are obtained. These equations can be rewritten 
as 
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The matrix M is composed by a diagonal part given by 
-^-Sij and by an off-diagonal part given by Ty. The 
diagonal part is only related with the material properties 
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of each atom, while the off-diagonal part only depends 
on the geometrical properties of the system. 

Once we solve the complex-linear equations shown in 
Eq. the dipole moment on each atom in the cluster 
can be determined. Then, we can find the extinction 
cross section of the whole nanoparticle, in terms of the 
dipole moments, as£& 

47rfc A 

Cext = -wr^Im(E* inc . Pl ), (4) 

i=i 

where * means complex conjugate. The extinction effi- 
ciency is defined as Q G xt = C cx t/^4, where A = Aird 2 , and 
d is the maximum distance between any two atoms in the 
nanocluster. Finally, the circular dichroism spectra is ob- 
tained by subtracting the averaged extinction efficiencies 
for left- and right-circularly polarized light, 

CD =< Qf xt > s - < Qcxt >s, (5) 

where < ■ • • >s indicates an average of the extinction 
efficiency, calculated over S different angular incident di- 
rections of the electric field^. 

In this work, we assume that the polarizability for each 
atom oti is the same for all the atoms in the cluster {pn 
— a), and is given by the well-known Clausius-Mossotti 
relation 

fR\ 3 (e-l) 

a ={V(7TTr < 6 > 

where e is the macroscopic dielectric function of the 
nanoparticle, and R is the ionic radii of the atom. The 
dielectric function is taken from the measurements by 
Johnson and Christyifi for bulk gold. This approximation 
has been used before providing good qualitative agree- 
ment between calculated optical properties and experi- 
mental data on gold and silver nanoparticlesii. On the 
other hand, since chirality is a geometrical property, the 
CD spectrum will be more sensitive to the second term 
inside the bracket in Eq. J3J , which depends on the clus- 
ter geometry, than to the first one that is related with the 
atomic polarizability. Therefore, in the present calcula- 
tions the simulated CD spectra will include the geomet- 
rical information about the chiral clusters, despite the 
polarizability is incorporated in an approximated form. 

The lowest-energy structures of the bare Au28 
and Au38, and thiol passivated Au28(SCH 3 )i 6 and 
Au38(SCH 3 ) 2 4 clusters had been obtained using density 
functional theory (DFT) within the generalized gradient 
approximation (GGA)i£. Both, the bare and thiol passi- 
vated most stable 28-atom gold clusters were found to be 
chiral, but a higher index of chirality was obtained for the 
passivated cluster. This lowest-energy passivated cluster 
was found upon relaxation of the bare chiral cluster, in- 
teracting with 16 SCH 3 molecules placed close to the 
three-atom hollow sites on the metal cluster surface^. 
The most stable thiol-passivated Au3s(SCH3)24 cluster 
was also found chirali*^, but it was obtained from the 



relaxation of the achiral AU38 cluster with Oh symmetry 
interacting with 24 SCH3 molecules placed on the three- 
atom hollow sites of the (111) faces of the bare cluster—. 
The AU38 isomer with Oh symmetry corresponds to a 
local minimum of the bare cluster according to DFT- 
GGA calculations^, while the lowest-energy structure of 
the bare AU38 was found with C s symmetry^, but upon 
passivation, it transforms into a higher-energy passivated 
isomer—. 

FIG. 1: Circular dichroism spectra of Au3g clusters, (a) Be- 
fore, with Oh symmetry, (dotted line) and after, with Ci sym- 
metry, (full line) passivation with 24 methylthiol molecules, 
(b) Lowest-energy isomer, with C s symmetry, (full line) and 
the same cluster geometry but taking out 10 atoms that are 
placed on the plane of symmetry (dotted line). 

In this work, we use the relaxed Cartesian coordinates 
of the clusters mentioned above to calculate their CD 
spectra. Figure I (a) shows the CD lineshape calculated 
using Eq. J5J) for the AU38 cluster before (dotted line) 
and after (full line) relaxation with the 24 methylthiol 
molecules. The zero value of the CD spectra before re- 
laxation indicates that the left- and right-circularly po- 
larized light is equally absorbed by the gold cluster since 
its geometrical structure corresponds to a truncated oc- 
tahedral structure with Oh symmetry. Since the geome- 
try of this cluster is highly symmetric, its CD spectrum 
should be zero in agreement with our result. On the other 
hand, the CD spectrum of the passivated cluster shows 
a non-zero lineshape due to the different absorption of 
the circularly polarized light. This behavior is expected 
since the passivated AU38 has been foun d chi ral using the 
geometrical approach described in Refs. Ill2l 

Figure 1 (b) shows an interesting result on the relation 
between chirality and the lineshape of the CD spectrum. 
The full line corresponds to the calculated CD spectrum 
of the lowest-energy structure of the bare Au38- The 
structure of this isomer has one plane of symmetry (C s 
symmetry) and therefore is achiralii£. However, as Fig. 
I (b) shows, the CD spectrum is non-zero. The reason of 
such behavior is coming from the atomic array existing 
on the plane of symmetry. This plane contains f atoms 
that are not in symmetric positions but forming a chi- 
ral two-dimensional array. This subset of atoms generate 
the non-zero signal of the CD spectrum. To confirm this 
effect, the CD lineshape of the 38- atom cluster without 
including the atoms on the plane of symmetry was calcu- 
lated (dashed line) . Since in this case there is a complete 
reflection symmetry in the cluster, its CD spectrum is 
negligible. 

FIG. 2: Circular dichroism spectra of AU28 clusters, (a) 
Before, with C\ symmetry, and (b) after, with C\ symmetry, 
passivation with 16 methylthiol molecules. 

Figures 2(a) and 2(b) display the CD spectra of the 
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Aii28 cluster before (bare) and after thiol relaxation (pas- 
sivated), respectively. Both clusters are chiral according 
to the geometrical analysis, based on the HCM, described 
in Refs. O]0 Therefore, it is expected that they have a 
non-zero signal in their CD spectra. However, Fig. 2 
also shows that indeed there exist notorious differences 
between the CD lineshapes of these clusters. These dif- 
ferences could account for the distinct indexes of chirality 
existing in these clusters as was calculated before using 
the HCMi£. In other words, Fig. 2 shows that clus- 
ters with different values of HCM can be distinguished 
through the lineshape of their CD spectra. 

To understand the origin of the differences in the CD 
spectra of clusters with different values of their indexes of 
chirality, the extinction efficiencies, using right-circularly 
polarized light, are shown in Figs. 3 (a) and (b) for the 
bare and passivated Au28 clusters, respectively. Since 
the CD spectrum is calculated from the difference of this 
quantity with the one obtained using left-circularly polar- 
ized light, it is interesting to analyze the most important 
features of its lineshape. From Figs. 3 (a) and (b) it is 
evident that the clusters with different HCM values also 
generate distinguishable extinction efficiencies. By ana- 
lyzing Eqs. 10-0), it is found that the cluster geometry, 
through the distribution of interatomic distances r^, is 
the main factor that generate the features (peaks and 
widths in the lineshape) of the extinction efficiency. 

FIG. 3: Extinction efficiency using right-circularly polarized 
light of Au28 clusters, (a) Before, with C\ symmetry, and 
(b) after, with Ci symmetry, passivation with 16 methylthiol 
molecules. The insets show the geometry of the bare and 
passivated clusters. 

The effect of the distribution of interatomic distances 
in the behavior of the extinction efficiency is more evi- 
dent by comparing this quantity for the 38-atom cluster, 
before and after relaxation with the thiol molecules. Fig- 
ure 4 shows the extinction efficiency using right-circularly 
polarized light for the bare AU38 isomer with Oh sym- 
metry (upper panel) and for the passivated chiral clus- 
ter (bottom panel). The insets in Fig. 4 display the 
distribution of interatomic distances of the correspond- 
ing clusters. From these figures it can be noticed that 
the extinction efficiency of the chiral passivated cluster 
(bottom panel) has more structure and shows a broader 



lineshape than the achiral cluster (top panel). This dif- 
ference can be correlated with the different behavior of 
the distribution of interatomic distances that indicate a 
discrete variation of distribution for the achiral Oh clus- 
ter, whereas the chiral one shows a much broader and 
continuous distribution of interatomic distances. 

The results presented in this work indicate that indeed 
chiral gold clusters with different indexes of chirality can 
be distinguished by the lineshape of their CD spectra. 
It was also found that within the dipole approximation, 
the origin of the distinct features of the CD lineshapes 
is related to the differences in the distribution of inter- 

FIG. 4: Extinction efficiency using right-circularly polarized 
light of AU38 clusters, (a) Before, with Oh symmetry, and 
(b) after, with C\ symmetry, passivation with 24 methylthiol 
molecules. The insets show the distribution of interatomic 
distances and the geometries of the bare and passivated clus- 
ters. 

atomic distances. Therefore, CD spectra are able to de- 
tect the geometrical features that characterize clusters 
with different indexes of chirality. These geometrical dif- 
ferences will also be important in quantum mechanical 
calculations of the CD lineshapes. It is encouraging that 
despite the approximations included in the calculation 
of the simulated CD spectra, their lineshapes show simi- 
lar features as those measured for glutathione-passivated 
gold nanoclusters^. In particular, there is a qualitative 
agreement between the calculated and experimental CD 
spectra in the energy range (2-4 eV), where the spectra 
shows the highest intensity . 

In summary, CD simulated spectra are appropriate for 
the characterization of clusters with different chirality 
and provide useful information for further studies using 
a quantum mechanical approach, and for experimental 
studies on the quantification of chirality in metal nan- 
oclusters. 
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